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We propose an experiment that would demonstrate nonlinear Josephson-type oscillations in the 
relative population of a driven, two-component Bose-Einstein condensate. An initial state is prepared 
in which two condensates exist in a magnetic trap, each in a different hyperfine state, where the 
initial populations and relative phase between condensates can be controlled within experimental 
uncertainty. A weak driving field is then applied, which couples the two internal states of the atom 
and consequently transfers atoms back and forth between condensates. We present a model of this 
system and investigate the effect of the mean field on the dynamical evolution. 



An interesting property of a weakly interacting Bose- 
Einstein condensate is that it can be ascribed an overall 
phase that can be measured relative to another conden- 
sate |I]-|e|. This is a quantum mechanical effect exhib- 
ited on a macroscopic scale. Some recent experiments on 
Bose-Einstein condensation (BEC) in dilute alkali vapors 
have investigated the relative phase of two overlapping 
condensates. For example, in the experiment reported in 
Ref. H , interference fringes in the density of two overlap- 
ping condensates were observed. More recently, the au- 
thors of Ref. [Q measured the relative phase of two con- 
densates in different hyperfine states using a technique 
based on Ramsey's method of separated oscillating fields 

A classic experiment that investigates the role of coher- 
ence on the dynamical evolution of two coupled macro- 
scopic quantum systems is the Josephson-junction exper- 
iment, in which a superconducting current of Cooper- 
pairs exhibits coherent oscillations Q. There have been 
several proposals suggesting an analogous experiment in 
the context of the more recent work being conducted on 
BEC of dilute atomic gases 

@|H- One can imagine 
preparing two initially isolated condensates in a double 
well potential and then lowering the central barrier to 
allow coupling between condensates to arise from tunnel- 
ing. One could then observe the time rate-of-change of 
the relative population, which is the analogous quantity 
to the current of Cooper-pairs in the usual Josephson- 
junction experiment. Further interesting effects could 
then also be studied, such as the nonlinear effect of the 
mean field on the system's behavior jf^|l(|. 

We propose a different experiment that would exhibit 
much the same physics as in the proposed double-well 
tunneling experiment but is based on the work done re- 
cently on two-component condensates, where two differ- 
ent hyperfine states can be populated and confined in the 
same trap J7|,p0|-p2[ . We restrict our attention to the sit- 



uation described in Refs. and (2l) where the authors 
trapped and cooled 87 Rb atoms in a magnetic trap below 
the critical point for BEC. The trapped atoms were ini- 
tially in the |/ = l,mf — —1} hyperfine state but after 
condensation the \f = 2, mj = 1} state could be pop- 
ulated through a two-photon transition. After this first 
7r/2-pulse, which lasts a fraction of the period of the trap, 
the relative motion of the two condensates oscillates and 
damps out to a stationary situation pl[ . After applying 
a second 7r/2-pulse, the authors observed a well defined 
relative phase that persists even beyond these damped 
oscillations J7|,pl|. 

In light of this observed "phase rigidity," we envision 
the following experiment. An initial stationary state is 
first prepared as described above by applying a short 
drive pulse that produces condensates in both 1, — 1} and 
2,1) states with known populations. When the tran- 
sient relative motion has damped out, the two conden- 
sates each sit in different shifted harmonic traps due to 
their different magnetic moments, with an overlap region 
that can be controlled experimentally. A weak driving 
field that couples these two internal states is then ap- 
plied so that the condensates are coupled in the overlap 
region. The time at which this sustained drive is turned 
on determines the initial relative phase accumulated by 
the condensates, which is measured relative to the ac- 
cumulated phase of the driving field 0. Oscillations in 
the relative population will occur, which depend on the 
initial relative phase and populations, and will exhibit 
nonlinear behavior due to the mean field. 

The true system described in Refs. and has ax- 
ial symmetry, with the two components separated along 
the z-axis axis. We simplify the problem by treating 
the system in only one dimension along the z-axis. We 
make another simplification by treating temperature zero 
and neglecting the fluctuations about the mean field that 
would give rise to damping. 
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In order to simplify our notation, we label the hyperfine 
states as |2) = |2, 1) and |1) = |1, —1). In the presence of 
a weak external magnetic field these hyperfine states are 
separated in frequency by ujq ■ The system is driven by a 
two-photon pulse, the strength of which we characterize 
by the two-photon Rabi frequency tt, which we take to 
be real- valued. We label the frequency of the two-photon 
drive by u)d, which can be varied to give a finite detuning 
5 = LJd — wq. In the following we have made the rotat- 
ing wave approximation by dropping the high-frequency 
terms in the atom-field interaction. Finally, we assume 
both states have a long lifetime compared to the period 
of the trap. 

After making the above approximations, we carried 
out standard mean-field theory on this coupled, two- 
component system to obtain the following equations for 
the time evolution of the condensates in the rotating 
frame: 

(i) 

The frequency of each trap is w z . We work in the "nat- 
ural" units of the problem, so that time is in units of 
1/uj z , energy is in terms of the trap level spacing ttu> z , 
and position is in units of the harmonic oscillator length 
z s ho = "v/Ti/wiRb oj z - The complex functions ipi(z,t) are 
the mean field amplitudes of each condensate, where 
i = 1,2. They are normalized to give the populations 
iVj(t), where the total number N — N 2 + N% is constant. 

The Hamiltonians appearing in Eq. ([!]) describe the 
free evolution and the mean field interaction 7J z MF 
for each component 

h™ f = xm 2 + K 3 m 2 , (2) 

where 71 = 1 and 72 = — 1, and z$ is the shift of each 
trap from the origin. The mean-field strength is char- 
acterized by Xij — fty/zsho, which depends on the scat- 
tering length dij of the collision. In general there will 
be three different values, one for each type of collision in 
this two-component gas: 022, an, 021- 

Before presenting results from numerical calculations 
of Eq. (0), it is useful to obtain two different forms of 
Eq. (0) that link this problem to two well known physical 
systems in the literature: the standard Rabi problem in 
quantum optics |p3| , and the Josephson-junction problem 
in condensed matter physics 0. 

We obtain the equations of motion for the popula- 
tions Ni and the coherences iVy = J dzijj*(z)ipj(z) from 
Eq. (|l|) by forming the appropriate products and inte- 
grating over space to yield 

N 2 = -iQ (Nn-Nia) 



N! = +in (N 2 i - n 12 ) 

N 21 = -iSN 21 +iA(t) - iQ{N 2 - Ni) , (3) 
where we define the time-dependent term A(i) as 

A(i) = -2z J dzz^(z)fa(z) 

+ Jdz (ff 2 MF - H™ F ) M(z)M*) ■ ( 4 ) 

The equations in Eq. (||) resemble the Bloch equations de- 
scribing an undamped, driven two-level atom |23) . How- 
ever, because the center-of-mass motion is correlated to 
the internal states of the atom, the extra term A(t) ap- 
pears, which includes the difference in external potentials 
between the two states. 

The first term in Eq. (^) arises from the difference in 
the shifted harmonic traps, which is just linear in z. This 
term acts as a time-dependent detuning. As population 
is transferred from one condensate to the other, the posi- 
tion of the overlap region changes due to the mean-field 
repulsion. This will cause the system to move in and 
out of resonance resulting in a suppression of the trans- 
fer of atoms. The second term comes from the difference 
in mean-field interactions and would vanish if all three 
scattering lengths were exactly degenerate, which can be 
seen from 

ffMF_ffMF = ( Aaa -A 21 )|V*| a 

- (An -A 21 )|</>i| 2 . (5) 

In order to make a link to the standard DC Joscphson 
effect, we must make some approximations in order to 
put Eq. (0) in a simpler form. For a very weak coupling 
(f2 << 1) and an initial state that is the self-consistent 
solution of the uncoupled system, we can make the ansatz 
ipi(z,t) = \jNi(t) e^*W $i(z). Here we put the explicit 
time dependence into the population Ni and the phase fa 
of each condensate while putting the spatial dependence 
into an adiabatic solution $i(z) to the undriven system 

(fl? + I?i" F )*i«=W*i(*), (6) 

where fii is the chemical potential for each condensate 
and the solutions &i(z) are taken to be real. The chemical 
potentials fa and functions vary slowly in time, 

being "slaved" by the populations. 

If we substitute this ansatz into Eq. (EI), we obtain the 
following equations of motion for the relative population 
V = (N 2 — Ni)/N and the relative phase tj> — (fa — fa) 

r) = -k(l -rf) 1/2 sin0 

j> = -K/z 2 - Ml) - S] + kr)(l - rfr^cos^ , (7) 

where k = 2 £1 J dz § 2 (z)$i(z) is proportional to the 
overlap of the condensates and so also varies slowly 
in time. These are non-linear versions of the usual 
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Josephson-junction equations M and are nearly identical 
in form to those obtained in Refs. |ll| and JlTj describing 
the double-well tunneling problem. The major difference 
is that in the double-well trap, the condensates are well 
separated, allowing the authors in Refs. and [jlTj to 
neglect the mean field in the interaction region of the 
barrier. In contrast, the interaction between condensates 
due to their significant overlap plays an important role 
in the evolution of the system described in this letter. In 
particular, it is this mutual interaction that causes the 
system to move out of resonance. 

We now show some results of calculations of both the 
exact solution given by Eq. ([!]) , and the approximate so- 
lution given by Eq. (Q) . The main purpose of the present 
calculations is to use realistic parameters to investigate 
the effect the mean field has on the evolution of the sys- 
tem. These parameters are listed in Table 1. 
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FIG. 1. This plot shows oscillations of the relative popula- 
tion. The dashed-dot line corresponds to the case where the 
mean-field interaction has been set to zero, the initial relative 
phase is 0(0) = n/2, and 5 — 0. The dashed line is for the 
same initial relative phase of n/2, but with the mean-field 
interaction turned on and 8 — —0.39. The solid line is for 
</>(0) = 37r/16 and S = —0.39. The three lines described were 
solutions of Eq. (Q) whereas the dotted line is a solution of 
Eq. (Q) with the same parameters as in the solid line. 

In Fig. 1 we show four curves that are described in 
the caption. As a point of reference, we plot the solution 
of Eq. (Q) with the mean-field terms set to zero. This is 
the standard Rabi solution, but here the Rabi frequency 
is given by cjr = 2 O J dz $2(z)&i(z), which includes 
the Frank-Condon-type overlap of the condensate wave- 
functions. However, when we turn on the mean-field in- 
teraction using the parameters given in Table 1, and set 
S = —0.39 so that the system is initially begin driven res- 
onantly, the amplitude is suppressed and the frequency 
has increased (the dashed line). As 0(0) is decreased, 
the amplitude increases, as illustrated by the solid curve 
where 0(0) = 3tt/16. Also, as 0(0) is decreased, the pres- 
ence of higher harmonics becomes stronger, as one can 



see in the shape of the solid line. 

We also plot the adiabatic solution given by Eq. for 
the case 0(0) = 37r/16 (the dotted line), where Eq. (||) 
is solved self-consistently in each time step. In this case 
the adiabatic solution agrees quite well with the exact 
solution (solid line) given by Eq. (g). The validity of the 
two-level, adiabatic solution depends on the structure of 
the evolving spectrum of this nonlinear system. In par- 
ticular, one must compare the time-rate-of-change of the 
Hamiltonian to the spacing between the instantaneous 
eigenmodes of the dressed basis Q. These quantities 
will depend on the size of the mean-field interaction, the 
strength of the coupling, given by 2f2 J dz ^2(z)^i{z), 
and also on the detuning S. It should also be noted that 
in the true system, the confining potential along the xy- 
plane is weaker than that along the z-axis, which should 
imply a more stringent criterion for adiabaticity than in 
the one-dimensional case considered. 
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FIG. 2. This plot shows the time-evolution of the density of 
each component. The x-axis is the position z, and the time t 
of each snapshot is shown. This case corresponds to the solid 
line plotted in Fig. 1. The detuning is <5 = —0.39, chosen so 
that the system is initially being driven on resonance. 

We plot a snapshot evolution of the densities of the 
condensates in Fig. 2 in order to show that the effect of 
the mean field is to push the system out of resonance as 
population is transferred between condensates. The case 
considered in Fig. 2 corresponds to the solid line in Fig. 
1. The detuning was chosen so as to compensate for the 
initial value of the term A(0) in Eq. (^), which represents 
the difference in external potentials, so that initially the 
system is being driven on resonance. However, as the sys- 
tem evolves, the first term in Eq. (g) gets larger since the 
|1) condensate is pushing the |2) condensate away from 
the center of the trap. This causes the region of over- 
lap to be displaced from the origin so that the system 
is no longer being driven on resonance. At t — 14.1 the 
region of overlap is centered at about z = 1.5, at which 
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time |A(£)/JV2i(£)| = 0.63, compared to the initial value 
|A(0)/JV 2 i(0)| = 0.39. This reduces the effectiveness of 
the drive and accounts for the suppression of the ampli- 
tude of oscillation in the relative population plotted in 
Fig. 1. 
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FIG. 3. This plot shows that the effect of the displacement 
zq on the system is to suppress the amplitude of oscillation. 
The four curves represent increasing values of zo: (solid), 
0.1 (dashed), 0.2 (dash-dotted), and 0.5 (dotted). In each 
case, the detuning S is chosen so that the system is initially 
being driven resonantly, and the initial phase is 0(0) = 7r/2. 

In Fig. 3 we vary the displacement z for four differ- 
ent values. The amplitude and period of oscillation both 
decrease as the displacement between the traps is made 
larger, while the overall shape does not vary much. In 
other words, the effective detuning caused by A(£) be- 
comes more pronounced for larger trap displacements. 
There is another effect of increasing zq: the region of 
overlap between the condensates decreases as the traps 
are pulled apart, which weakens the coupling. One might 
expect the period of oscillation to increase in this case, 
however, the plot displays a decreasing period as zq is 
increased. This indicates that the time-dependent de- 
tuning due to A(£) has a stronger effect. 

In this letter we have described a physical system based 
on the experiments reported in Refs. |20|j2l| j7|] that would 
exhibit non-linear Josephson-like oscillations in the rela- 
tive population between a driven two-component conden- 
sate. From our calculations we have observed the effect 
of the mean field, which acts to suppress the transfer of 
atoms between overlapping condensates in separated har- 
monic traps and gives rise to nonlinear oscillations in the 
relative population. We have also found that, for equal 
populations initially 77(0) = 0, the spectrum of Fourier 
components comprising the oscillation depends on the 
initial relative phase </>(0). 
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cussions. We also thank the National Science Foundation 
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TABLE I. This is a table showing the values used for the 
various physical parameters appearing in our calculations. 
The scattering lengths are taken from Ref. [21]. Since we are 
treating the system in ID, N is not the actual population, 
but was chosen to produce a realistic mean-field interaction 
for 5 x 10 5 atoms [7]. 
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